Anthropic predictions for vacuum energy and neutrino masses in the light of WMAP-3 
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Anthropic probability distributions for the cosmological constant A and for the sum of neutrino 
masses m„ are updated using the WMAP-3 data release. The new distribution for A is in a better 
agreement with observation than the earlier one. The typicality of the observed value, defined as the 
combined probability of all values less likely than the observed, is no less than 22%. We discuss the 
dependence of our results on the simplifying assumptions used in deriving the distribution for A and 

• show that the agreement of the anthropic prediction with the data is rather robust. The distribution 

for m„ is peaked at m v ~ 1 eV, suggesting degenerate neutrino masses, but a hierarchical mass 

, pattern is still marginally allowed at a 2o~ level. 
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<D , I. INTRODUCTION 

Q. 

The parameters we call constants of Nature may in fact be stochastic variables taking different values in different 
parts of the universe. The observed values are then determined by chance and by anthropic selection. Recent 
developments in string theory [1, 2], combined with the ideas of eternal inflation [3, 4], have led to the "landscape" 
paradigm, which provides a theoretical basis for this picture. And the successful prediction of a nonzero cosmological 
constant A [5-10] may be our first evidence for the existence of a huge multiverse beyond our horizon. (For an up to 
date discussion of these ideas, see [11].) 

Anthropic predictions are necessarily statistical in nature. Even though the observed value of A is within the 
expected range, there is a lingering doubt that this may be no more than a lucky guess. To further test the theory, 
we should try to extract predictions for variables other than A. Some attempts in this direction have already been 
made. Predictions for the neutrino masses m v have been discussed in [12, 13], and a postdiction for the (axionic) 
dark matter mass per baryon, £, has been given in [14]. 

A reliable calculation of anthropic probability distributions is a very challenging task. All calculations performed so 
O far have relied on a number of simplifying assumptions. Moreover, the form of the distributions for A and m„ depends 
JL \ on the values of the cosmological parameters: the Hubble parameter h, the tilt of the perturbation spectrum n s , and 
; I ■ particularly on as - the amplitude of linear density perturbations on the scale of 8ft." 1 Mpc. Hence, the agreement 
(or disagreement) of the predictions with the data should be regarded as tentative. With improved understanding or 
Ci • improved measurements, it may either get better or worse. 

In Fig. 1 we show the distribution P(A) circa 2005. Its calculation has been discussed in [7-10, 15-17]. The 
• i— i , parameter values used in the figure were based primarily on the first year of WMAP data. The agreement between 
the theory and observation can be characterized by the "typicality" of the observed value A*, T(A») [18]. (The 
definition of typicality will be discussed in Section II.) For the distribution in Fig.l, T(A«) » 10%. For a statistical 
. - - i model of this sort, this can be regarded a good agreement. [48] 

With the release of the 3-year WMAP data, some of the parameters (such as n s and h) are now known with a 
better accuracy. Among implications was a decrease in the value[49] of as by more than 10%, from 0.9 ± 0.1 [19] to 
0.76 ± 0.05[20]. Combining other datasets with the WMAP data can change the preferred values of the parameters. 
For example, adding all available information, including the SDSS Lyman-a forest data [21], was reported to increase 
the value of a 8 to 0.85 ± 0.02 [22]. On the other hand, combining the WMAP data with the distribution of SDSS 
luminous red galaxies yields as = 0.76 ± 0.03 [23] in good agreement with results of the WMAP team. Here, as in 
our original analysis, we will make use of the WMAP data alone, which is most likely to be free of systematics and, 
for the flat power law ACDM model, provides sufficiently accurate measurements of the spectral index and other 
cosmological parameters that determine erg. 

The main goal of the present paper is to update the distributions for A and m„ using the WMAP-3 data. We 
shall see that both distributions are shifted toward smaller values of A and m„, respectively. The agreement between 
theory and observations is improved as a result. We shall also discuss the dependence of the results on the simplifying 
assumptions made in the calculation of probabilities. 

In the next section, we discuss the concept of typicality which will be used as a quantitative measure of the agreement 
of theoretical distributions with the data. The case where A is the only variable is considered in Section III. The 
effects of other variable parameters, such as the amplitude of the primordial density perturbations Q and the dark 
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FIG. 1: The probability distribution for A based on the flat ACDM best fit parameters from WMAP-1. A* is the observed value whose 
la uncertainty is represented by the blue shaded vertical stripe. 

matter mass per photon £, and the dependence on the assumed value of the galactic mass Mq, are discussed in IV. 
Variable m„ and joint variation of m v and A is analyzed in Section V. We finish with some concluding remarks in 
Section VI. 



II. TYPICALITY 

Consider a variable X defined in the interval a < X < b and characterized by a probability distribution P(X). The 
value X — X* has typicality T if T% of the distribution are "less probable" than X*. Page [18] assumes that the 
least probable values are at the tails of the distribution (see Fig. 2). The motivation is that such values are fine-tuned 
to be near the special points a and b, which are the endpoints of the interval. 

Suppose for definiteness that X* is closer to a than to b in terms of the measure P(X), so that 

p* = £' P(x)dx < 1 (l) 

Let us also define Xj, which is equally close to b, 

[ P(X)dX = P*. (2) 

The values in the intervals a < X < X* and Xj < X < b are regarded as less likely than V*. The typicality is then 
given by 

T(X0 = 2P*. (3) 

In the case of the cosmological constant, apart from the vicinity of the endpoints at A = ±oo, there is an additional 
range of unlikely values near zero. A — is a special value, which defines the boundary between eternal expansion 
and recollapse. Having A accidentally fine-tuned very close to this point is unlikely, and we want this reflected in the 
typicality. We shall define P* as the "distance", in terms of the measure P(A), from the observed value A* to the 
nearest special point (A = 0, ±oo). It will also be convenient to define the probabilities P± for positive and negative 
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FIG. 2: Values "less probable" than X* are located in the two shaded regions, which have equal measure in the distribution P(X). 

A, 

P+ = / P(A)dA, (4) 
Jo 



-f 



P(A)dA, (5) 



P + + P_ = 1 . A straightforward generalization of Page's definition is to define typicality as the measure of all points 
whose distance to the nearest special point is less than P*. 

It is clear from Fig. I that the special point closest to the observed value A* is A = 0. Hence, 



P* = / * P(A)dA < P+/2. (6) 
Jo 

If in addition P* < P_/2, then there are four intervals of A less likely than A*, each having measure P*, and thus 

T(A,) = 4P,. (7) 

Alternatively, if P* > P_/2, we have[50] 

T(A») = P- +2P*. (8) 

Most of the calculations of the distribution P(A) have been performed for A > 0. For negative A, the universe 
recollapses at t\ — (67rG|A|) -1 / 2 . To find the distribution in this case, we need to know the probability for a 
civilization to evolve during this time interval, which is of course very uncertain. It is very likely, however, that 
positive A are more probable than negative A [24] , 

P+>P- (9) 



Note, for example, that for A = — 20A* the universe would recollapse in 3 billion years. Star formation in such a 
universe will not begin before the onset of the recollapse, and the remaining 1.5 billion years are hardly sufficient 
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to disperse heavy elements in supernova explosions, form planetary systems and evolve intelligent life. On the other 
hand, the 2a range for positive A extends to 80A„. The inequality (9) can be used to derive a lower bound on T(A*). 
Let 7+(A») be the typicality calculated in the interval < A < oo with P(A) normalized in that interval, 

T+(A») = 2/yP+. (10) 

Now, using Eq. (9) it is easily verified that 

T(A„) > T + (A»). (11) 

In the rest of the paper we shall consider only positive A and the corresponding typicality 7+. According to (11), 
it can be regarded as a lower bound for the full typicality T. 



III. THE PROBABILITY DISTRIBUTION FOR A 



The method used in this paper follows closely that in [10, 12, 13, 16]. Here we only reproduce the essential details 
and refer the reader to [13] for a comprehensive description. Consider a model in which A is allowed to vary from one 
part of the Universe to another. We define the distribution V(A)dA as the probability for a randomly picked observer 
to measure A in the interval dA. This distribution can be represented as 

P(A) = P pnor (A)f selec (A), (12) 

where the prior probability P pr i or (A) is determined by the geography of the landscape and by the dynamics of eternal 
inflation, and f se iec(A) accounts for anthropic selection effects. 

The standard argument [9, 15] suggests that the prior probability is well approximated by 

^^'pviov i ™° const . ^13^ 

because the anthropic range where f se iec(A) is appreciably different from zero is much narrower than the full range 
of A. We emphasize that this is just a heuristic argument. The conditions for its validity have been discussed both 
in scalar field models, where A is a continuous parameter [16, 25-27], and in "discretuum" models with a large set of 
metastable vacua [28, 29]. Future work will show whether or not these conditions are satisfied in the string theory 
landscape. Here we shall assume that Eq.(13) is valid. 

It has been the standard practice to identify the selection factor in (12) with the asymptotic fraction of baryonic 
matter, F(M > Mq,A), which clusters into objects of mass greater than the characteristic galactic mass Mq <~ 
1O 12 M : 

/ seiec (A)cxF(M>M G ,A). (14) 

The idea here is that there is a certain average number of stars per unit baryonic mass and a certain number of 
observers per star, and that these numbers are not strongly affected by the value of A. In the present Section we shall 
adopt this standard approach; its validity will be analyzed in Section III. 

The fraction of collapsed matter F(M > Mq,A) can be approximated using the Press-Schechter (PS) formalism 
[30]. This leads to 



F(M > Mq, A) oc crfc 

where 



X 



V2 



(15) 



X = 6 /a oo (A) , (16) 

S c rts 1.63 and (Too is the variance of the Gaussian density fluctuation field in the asymptotic future on the galactic 
scale Mq. The latter quantity can be written as [13] 

'-< A >-*»5(E^)- (17) 

Here, a* G is the current density contrast on the galactic scale Mq inferred from the large-scale CMB data, x % is our 
local value of Q^/Q M , A* is the local value of vacuum energy density, D(A*, x») is the local linear growth factor and 
D°°(A) = D(A,x = oo). 
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N-body simulations indicate that the PS model overestimates the abandance of "typical" halos, while underesti- 
mating that of more massive structures [31]. The Sheth-Tormcn model [32] was shown to fit the simulations better. 
We have checked that replacing the PS formula with that of [32] does not significantly affect our results and use the 
PS model through out this work. 

To evaluate u* G we first find the length scale R(M G ) corresponding to the mass scale M G using 



R(M G ) 



/3M G \ 1/3 



0.951 



M 



1/3 



Mr 



10 12 M P 



1/3 



h^Mpc , 



(18) 



where p ~ 1-88 • 10 26 VL M h 2 kg/m 3 is the mean cosmic density of nonrelativistic matter. For h = 0.73 and 
CIm = 0.234 this gives 



R(M G ) w 1.4 h _1 Mpc 



10 12 M m 



1/3 



(19) 



The corresponding linearized density contrast a* G found using WMAP-3's best fit power law model [20] for M = 
1O 12 M is 



a* G = a[R(M G ) 



1 

2^2 



P{k)W 2 {kR)k 2 dk 



1/2 



1.87 ±0.13 . 



(20) 



In the above, P(k) is today's linear matter power spectrum and W(kR) is the window function. We work with the 
"top-hat" form for W(kR), also used by the WMAP team. This value is smaller than a* G w 2.41 ± 0.26 inferred from 
the WMAP-1 data and used in [12, 13]. 

The growth factor D(A, x) in a universe containing only A and pressureless matter can be written as [10, 33] 



D(A,x) = l + -x eq -^ 3 G(x) 



(21) 



where 
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1/2 r dw 
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X 


J itfVe (l + w )3/2 



(22) 



Using \x eq 1 / 3 G(x) » 1, the variable X defined in Eq. (16) can be written as 



x = -\ A 



1/3 



A*/ G(oo) 



(23) 



The variance a* G is proportional to the amplitude of primordial fluctuations Q. The dependence on Q can be 
explicitly introduced by writing 



* Q 



(24) 



where Q* is the observed value of Q. We follow the WMAP-3 team's conventions [20] and define Q = Ak(k 
0.002h/Mpc) with the best fit value of Q* = 4.9 x 10" 5 . Then X can be written as 



X = A 



1/3 



Q ' 



with 



A w 3.2 x 10" 5 



10 12 M n 



o.i 



(25) 



(26) 



where the mass dependence of A comes from that of a G in (20) and the approximation is valid for values of M G within 
a few orders of magnitude of 10 12 M©. 
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FIG. 3: The probability distribution for A obtained as in Fig. 1 but using parameters derived from the WMAP-3 data. The observed 
value, A*, has a typicality of 22%. 



In summary, the distribution for A is 

" X 



P(A)dA oc erfc 



V2. 



dk, (27) 



where X is given by Eq. (25) with A from Eq. (26). 

The degree of clustering that we observe today is often parameterized by the variance of the linear density contrast 
field on the 8/i _1 Mpc scale, a%. The latest observations, as mentioned in the introduction, indicate a lower value of 
as, and correspondingly a lower value of oq- A decrease in ac at a fixed A leads to a suppression of galaxy formation, 
except for very small values of A, such that A-domination occurs when clustering on the galactic scale is essentially 
complete. As a result the peak of the probability distribution is shifted towards lower values, leading to a better 
agreement with the observed A. 

The updated distribution is shown in Figure 3 and corresponds to a typicality of 7+ = 22% for our local value of A. 

In [24] Weinberg reports a probability of 15.6% for the observed value of A, implying a typicality 7+ = 31%. (His 
probability includes only values between and A*; it is half of what we call 7+.) This result is based on the WMAP-1 
data and is higher because of the use of the Gaussian window function with Rq = 2 Mpc. This corresponds (for the 
same mass) to an effective "top-hat" radius of R(Mg) ~ 2.3 Mpc/h (as opposed to 1.4 Mpc/h used in this paper). 
For the WMAP-3 parameters Weinberg's calculation would ycild a typicality around 50%. 



IV. DEPENDENCE ON SIMPLIFYING ASSUMPTIONS 



The calculation of P(A) that led to the distribution in Fig. 3 relied on a number of simplifying assumptions. In 
particular, we assumed (i) that A is the only variable parameter and (ii) that the number of observers is proportional 
to the mass fraction in giant galaxies of mass greater than Mq = 10 12 M Q . We shall now discuss how sensitive our 
conclusions are to these assumptions and what modifications we can expect in more realistic models. 



A. Other variable parameters 



A criticism that has often been raised against the anthropic prediction of A is that the agreement of the prediction 
with the data is destroyed if one allows parameters other than A to vary. For example, if the amplitude of the 
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FIG. 4: The probability distribution for A assuming three different mass thresholds for halos that can become habitable galaxies. From 
left to right, the curves correspond to thresholds of 10 12 Mq, 10 11 Mq and 10 10 Mq . The typicalities of the observed value, A*, are 22%, 
11% and 7%, respectively. 

primordial density perturbations Q is also variable, then galaxies will form earlier in regions with larger Q, and A 
can take larger values without interfering with galaxy formation [34]. However, it has been shown in [35] that the 
probability distribution in this case factorizes as 

P(Q,A)dQdA = Pp„ or (Q) erfc 

Here the factors Q 3 and X 2 come from the transformation Jacobian and we have used the expression for X in Eq. (25). 

The value of Q is determined by the shape of the inflaton potential, and the prior distribution for Q is highly model- 
dependent [36-38]. So it is very fortunate that the Q-dependent part of the distribution has completely factored out, 
and we are left with the distribution for X, 

d(X 3 ). (29) 

For a fixed Q, X 3 cx A with a constant coefficient, and we recover Eq. (27). For a variable Q, Eq. (29) shows that the 
distribution has exactly the same form, except now it should be regarded as a distribution for the variable X 3 cx A/Q 3 . 
The typicality of the observed value of X is therefore the same as what we found for A* in the preceding section, 
T(xJ) = 22%. 

This conclusion can be extended to include a variable dark matter mass per photon, £. Using the analysis in [14], it 
can be shown that, to a good approximation, the distribution in this case can still be factorized into a part depending 
on Q and £ and a part depending on X 3 cx A/Q 3 £ 4 . The predicted range of X is unaffected by variation of Q and £ 
and is independent of their prior distributions. 

B. Dependence on the characteristic mass scale Mg 

Next we investigate the dependence of P(A) on the galaxy mass cutoff Mg. It has been argued by Loeb [39] that 
the anthropic explanation of A would be significantly weakened if a large number of habitable planets were found to 
exist in smaller galaxies, such as the dwarf descendants of galaxies formed at z ~ 10. It would mean that planet-based 
observers could be abundant in our Universe even if the cosmological constant was some orders of magnitude larger 
than observed. 




X 

.71 



dQdA oc Q 3 P P rior{Q)dQ erfc 



X 

.71 



X dX. 



(28) 



P(X)dX cx erfc -= 
_v2_ 
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FIG. 5: The 68 and 95% joint probability contours for A and m„ obtained using the WMAP-3 parameters. 



One motivation for introducing the cutoff mass Mq is that, in the hierarchical structure formation scenario, smaller 
galaxies typically form at earlier times and have a higher density of matter. This may increase the danger of nearby 
supernova explosions and the rate of near encounters with stars, large molecular clouds, or clumps of dark mater 
[14, 16, 40]. Gravitational perturbations of planetary systems in such encounters could send a rain of comets from the 
Oort-type clouds toward the inner planets, causing mass extinctions. (Encounters close enough to disrupt planetary 
orbits are much less likely.) With the present rate of mass extinctions (about once in 10 s yrs), the intervals between 
these events are already comparable to the time it took humans to evolve. Hence, a substantial increase of the rate 
may result in a suppression of the number of observers. 

Another consideration is that the properties of galaxies, as a function of their mass, are observed to undergo a sharp 
transition at the halo mass of about Mh ~ 2 x 10 n M© (stellar mass ~ 3 x 1O 1O M0) [41]. For larger masses, baryons 
are efficiently converted into stars, while for smaller masses the efficiency is lower and decreases as M gets smaller. 
The metallicity is also observed to drop with mass [42]: dwarf galaxies appear to be deficient in heavy elements, which 
are necessary for the formation of planets. These traits are attributed to supernova feedback: supernova explosions 
either expel a fair fraction of gas from low-mass galaxies or inhibit star formation. (For a discussion and references, 
see [43]). Low mass fraction in stars and low metallicity indicate that dwarf galaxies are less likely to harbor observers. 

The effective cutoff mass Mq cannot be determined without a quantitative study of these effects. We shall not 
attempt this here. In order to quantify the dependence of our predictions on the value of Mg, we repeated the analysis 
in Section III using two lower values, namely, 10 11 Mq and 10 10 Mq. The resulting distributions are shown in Fig. 4. 

As expected, a decrease of Mq shifts the distribution towards larger values of A. However, the effect is not 
dramatic. The observed value remains in the 95% range of the distribution. Its typicality 7+ is reduced from 22% for 
Mq = 10 12 M Q to 11% and 7% (for Mq = 10 n Af o and Mq = 10 10 M Q , respectively). 



We now consider the case when both the sum of the neutrino masses, m v = m\ + m,2 + fn^, and the vacuum energy 
density, A, vary from one part of the Universe to another. The possibility that the smallness of the neutrino masses 
may be due to anthropic selection has been first suggested in [12]. The mass-squared differences of the three neutrino 
species are now known with an accuracy of about 10%, 



V. 



VARIABLE m„ AND JOINT VARIATION OF m„ AND A 



Am 2 ,. 



?,2 



2.6 x 10~ 3 eV 2 , 



(30) 



Amlj « 8 x 10~ 5 eV 2 . 



(31) 
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FIG. 6: The probability distribution for the sum of the neutrino masses, m„, obtained using the WMAP-3 parameters. 



The absolute values of the neutrino masses are currently unknown. This makes them an ideal target for anthropic 
predictions. 

It is usually assumed that the neutrino mass eigenvalues form a "normal" hierarchy, i.e., mi « m2 < «i3. Then 
it follows from Eqs. (30), (31) that to„ <~ 0.06 eV. The opposite regime is when the three masses are degenerate, 
mi ss rri2 ~ 7713. An intermediate case is that of inverted hierarchy, 771,3 ^ m i ~ mi, m which case m„ « 0.11 eV. The 
astrophysical upper bound from the WMAP-3 data is m u < 2 eV. Stronger bounds have been quoted from combined 
astrophysical data sets, but their validity is not certain (for a recent review see [44]). 

Assuming that m„ and A vary independently, their joint prior probability is a product of the individual priors. We 
have already chosen a flat prior for A in Section III. The prior probability for m v was discussed in [12, 13] with the 
conclusion that V pr ior{m v ) — const is the best motivated choice. As in Section III, we approximate the anthropic 
selection factor as the asymptotic mass fraction in galaxies of mass M > Mq, 



f se iec(m v ,A)ocF(M>M G ,m l/ ,A) oc erfc 

where 



X 



V2 



(32) 



y ^ c D\f v *, A*, a;*) ,„„> 

X 'a(M) £><»(/„, A) • (33) 

The notation above is similar to that of Eq. (17) of Section III, except for the hats indicating the quantities evaluated 
for our local region assuming massless neutrinos. We refer the reader to [13] for a detailed derivation of the above 
expressions. The updated joint probability distribution for A and 777„ is shown in Figure 5. 

It is of interest to also evaluate the anthropic prediction for m„ under the assumption of fixed A = A*. This 
distribution is shown in Figure 6. The difference, compared to the results in [12], is that the mean value of m„ has 
decreased from 2.1eV to 1.5eV, and the lower boundary of the 2a region from 0.07 eV to 0.05eV. 

Both the single-variable and combined distributions are peaked at a value m„ <~ 1 eV, corresponding to degenerate 
neutrino masses. The la range of the distribution, 0.3 eV < m v < 3 eV, also predicts degeneracy. The prediction, 
however, is not very strong: the 2cr range is (marginally) consistent with the hierarchical mass pattern. We note that 
for 777^ <~ 1 eV the observed value of A is well within the la range of the distribution. 
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VI. DISCUSSION 



The observed value of the cosmological constant, A*, is mysterious for two reasons. First, all particle physics 
estimates yield enormous values, many orders of magnitude greater than A,. Second, A, is comparable to the average 
matter density at the present epoch - a coincidence calling for an explanation. Both of these facts are accounted for 
in the multiverse model, where A is a variable changing from one place in the universe to another. Despite many 
attempts, there are no plausible alternative explanations. A satisfactory agreement between the observed A* and the 
anthropic prediction is therefore hard to dismiss. It may be our first evidence for the existence of the multiverse. 

In this paper we have updated the anthropic probability distribution for A using the WMAP-3 data release. This 
resulted in an improved agreement between the theory and the data. The typicality of the observed value A* , defined 
as the probability of all values "less likely" than A*, is 7+(A„) = 22%. (As we explained in Sec. II, this estimate 
should be regarded as a lower bound, since our analysis did not include negative A.) 

We analyzed the dependence of the successful prediction for A on the simplifying assumptions used to derive the 
distribution P(A). Our conclusion is that the prediction is rather robust. It is not drastically modified by variation 
of the cutoff galactic mass Mq and remains essentially unaffected by inclusion of other anthropic variables, such as 
the amplitude of primordial density perturbations Q, or mass of dark matter per baryon £. Variation of the sum of 
the neutrino masses vn v docs have some effect on P(A). As m„ is increased from its lowest value, the agreement with 
the data improves, and for m v ~ 1 eV the observed value A* is well within the la range of the distribution. 

We have also updated the anthropic prediction for m v . With a flat prior, the distribution is peaked at m v ~ 1 eV, 
suggesting degenerate neutrino masses. A hierarchical neutrino mass pattern is marginally acceptable at the 2a level. 

Unlike the case of the cosmological constant, there is a viable alternative explanation for the smallness of the 
neutrino masses. It is the see-saw mechanism. We note that a recent study failed to find models exhibiting this 
mechanism in a wide class of string theory inspired models [45]. We note finally that a value of m v ~ 1 eV is 
suggested by the Heidelberg-Moscow double beta decay experiment [46]. This claim, however, remains controversial 
[47]. 
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